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ABSTRACT 

The environment surrounding Wolf-Rayet (W-R) star HD 211853 is studied in 
molecular, infrared, as well as radio, and Hi emission. The molecular ring consists of 
well-separated cores, which have a volume density of 10^ cm"^ and kinematic temper- 
ature ~20 K. Most of the cores are under gravitational collapse due to external pressure 
from the surrounding ionized gas. From the spectral energy distribution modeling to- 
ward the young stellar objects, the sequential star formation is revealed on a large 
scale in space spreading from the W-R star to the molecular ring. A small-scale se- 
quential star formation is revealed toward core "A", which harbors a very young star 
cluster. Triggered star formations are thus suggested. The presence of the photodisso- 
ciation region, the fragmentation of the molecular ring, the collapse of the cores, and 
the large-scale sequential star formation indicate that the collect and collapse process 
functions in this region. The star-forming activities in core "A" seem to be affected by 
the radiation-driven implosion process. 

Subject headings: ISM: bubbles C ISM: molecules C stars: formation C stars: Wolf-Rayet 



1. Introduction 

The feedback of massive stars is expected to strongly influence their surrounding 
interstellarmedium (ISM) through stellar radiation, energetic wind, and the ejection of the 
elemental composition. The shock front (SF) that emerges as bubbles around massive stars expand 
can compress the ISM and trigger star formation in very dense layers (Chen et al. 2007). 

Two models have been proposed to explain the effect of massive stars on star formation: 
collect and collapse (Elmegreen & Lada 1977) and radiation-driven implosion (RDI; Lefloch & 
Lazareff" 1994; Ogura 2010). The former takes place in a larger spatial size (~10 pc) with a longer 
timescale (a few Myr), whereas the latter takes place in a smaller spatial size (~1 pc) with a shorter 
timescale (0.5 Myr; Ogura 2010). In the "collect and collapse" scenario, the ionization front 
(IF) generated from the Hii region gathers molecular gases and collects a dense shell between 
the IF and the SF. The gas and dust in the collected layer collapse to form stars when they reach 
the critical density (Chen et al. 2007; Ogura 2010). This process can selfpropagate and lead to 
sequential star formation (Elmegreen & Lada 1977). Star formation induced through the "collect 
and collapse" process has been revealed in the borders of several Hii regions (Deharveng et al. 

2003, 2005, 2008; Zavagno et al. 2006, 2007; Pomar'es et al. 2009; PetrieUa et al. 2010;Brand 
et al. 2011). The so-called RDI process has been proposed to account for the evolution of 
bright-rimmed clouds whose surfaces are ionized by the UV photons from massive stars and 
trigger the star formation therein (Bertoldi 1989; Bertoldi & McKee 1990; Hester & Desch 2005; 
Miao et al. 2009; Gritschneder et al. 2009; Bisbas et al. 2011; Haworth & Harries. 2012). The IFs 
generate a shock that compresses into the cloud and causes the dense clumps to collapse (Chen et 
al. 2007). Candidates for star formation by RDI have been found in many bright rimmed clouds 
associated with IRAS point sources (Sugitani et al. 1991; Sugitani & Ogura 1994; Urquhart et al. 

2004, 2006, 2007; Morgan et al. 2010). 

The Wolf.Rayet (W-R) stars, initially more massive than ~40Mo, have typical mass-loss 



rates of 10"^ M© yr"' and strong stellar winds faster than 1000 km s ^ They can blow stellar 
wind bubbles with radii of ~10 to ~100 pc and expanding velocities of several km s"' (Marston 
1996;Martin-Pintado et al. 1999). The expanding bubbles can push away the molecular gas and 
create molecular layers that fragment into condensations well separated in due time. The "collect 
and collapse" process is expected to happen in the shells swept up by such expanding stellar wind 
bubbles (Whitworth et al. 1994). However, the samples of triggered star formation surrounding 
W-R stars are rare and few previous works have focused on star-forming activities in such regions. 
In this paper, we provide a sample region where triggered star formation is witnessed. 

HD 21 1853 is located at the center of an optically bright region named shell B in the Hii 
region Sh2-132 (Sharpless 1959; Harten et al. 1978). The Hii region linked to shell B is excited 
by a massive 0-type star BD+55.2722 and the W-R starHD21 1853 (Vasquez et al. 2010). While it 
is the main ionizing source of the ring nebula (Cappa et al. 2008), the W-R star HD 21 1853 (WR 
153ab) is a binary system. The spectral types of the two W-R components are WN6 and WCE-i-06 
(Nishimaki et al. 2008). The velocity of its stellar wind is as high as 1500.1800 km s.l, leading to 
a high mass-loss rate of 10"*^ Mq yr"^ (Nishimaki et al. 2008). A ring nebula associated with HD 
21 1853 is easily identified in radio emission, in polycyclic aromatic hydrocarbon (R^.H) emission, 
and in molecular emission (Vasquez et al. 2010; Cappa et al. 2008). Figure 1 presents a composite 
color image of this region. Blue corresponds to the DSS-R optical emission, green to the radio 
emission detected in the Canadian Galactic Plane Survey (CGPS; Taylor et al. 2003) at 1420 
MHz, and red to 8.3 yum emission detected at the Midcourse Space Experiment (MSX) A band. 

It can be seen that the radio emission is embraced by the 8.3 /urn emission, which forms a 
ring-like structure. A photodissociated region (PDR) has been identified, and possible triggered 
star formation in this region is suggested (Vasquez et al. 2010). As shown in Figure 2 of Saurin 
et al. (2010), three star clusters, "Teutsch 127," "SBB 1," and "SBB 2," have been found in 
this region; the hierarchical structure and age distribution of the star clusters show evidence of 
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triggered star formation (Saurin et al. 2010). However, studies of the star-forming activities in 
this region are far from mature. In this paper, through detailed analysis of molecular emission 
and characteristics of young stellar objects (YSOs), we have revealed the picture of triggered star 
formation in this region. 

2. Observations and Database 

The observations of ^^CO (1-0), "CO (1-0), and C^^O (1-0) were carried out with the 
Purple Mountain Observatory (PMO) 13.7 m radio telescope in 2011 January. The new 9-beam 
array receiver system in single-sideband (SSB) mode was used at the front end. The array 
pixels are sideband separation SIS mixers yielding 18 independent IF outputs. The prop erties of 



Shan et al. 



( 2009h . i^CO 



the Solid-State Imaging Spectrometer (SIS) mixers are introduced in 
(1-0) was observed at the upper sideband, while "CO (1-0) and C^^O (1-0) were observed 
simultaneously at the lower sideband. FFTS spectrometers were used as back ends, which have 
a total bandwidth of 1 GHz and 16384 channels, corresponding to a velocity resolution of 0.16 
km s' for "CO (1-0) and 0.17 km s"' for "CO (1-0) and C^^O (1-0). The half-power beam 
width was 56", and the main beam efficiency was ~ 0.5 during our observation period. The 
pointing accuracy of the telescope was better than 4". The typical system temperature (T^^^) in 
SSB mode was around 110 K and varied about 10% for each beam. The On-The-Fly (OTF) 
observing mode was applied. The antenna continuously scanned a region of 22' x 22' centered 
on the Wolf-Rayet star HD 211853 with a scan speed of 20" s ^ The OTF data were then 
converted to three-dimensional cube data with a grid spacing of 30". Because the edges of the 
OTF maps are very noisy, only the central 14' x 14' region was selected to be analyzed. The 
typical rms noise level was 0.2 K in T;; for "CO (1-0), and 0.1 K for "CO (1-0) and C'*^0 (1-0). 
] ^or the data analysis, the G ILDAS software package including CLASS and GREG was employed 



( Guilloteau & Lucas 



20001). 
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Neutral hydrogen (HI) 21 cm line data were from CGPS with intensities shown as brightness 
temperature T,,,/^. The synthesized beam and rms bright temperature in a 0.82 km s"^ channel 
were ~ 1' and 3 K, respectively. R adio continuum data and images at 1420 MHz were extracted 



from the CGPS (ITaylor et al 



Condon et al. 



2003h and the NRAO Very Large Array (VLA) Sky Survey (NVSS, 



(119981)'). The positio ns and fluxes of the NVSS sources were obtained from the 



NVSS catalog (.Condon et al. 



119981) . The CGPS 1420 MHz data had a synthesized beam of ~ 50" 
and a noise level of ~ 3 mJy beam"' . The NVSS image had a spatial resolution of 45" and an rms 
brightness fluctuation of ~ 0. 45 mJy beam"' (S toke I). MSX A-band data were extracted from the 



MSX Galactic Plane Survey (iPrice et al. 



20011), which has a spatial resolution o f ~ 18". 3 



IRAC data and Two Micron All Sky Survey (2MASS) data were obtained from 



Oiu et al. 



Spitze r 



(l2008h . 



both of which have good spatial resolutions of ~ 2". We also obtained SCUBA 450 jum and 850 



yum continuum data from the SCUBA Legacy Catalogues (IDi Francesco et al. 



20081). At 850 and 



450 yum, the resolutions of SCUBA data are of ~ 14" and ~ 9", respectively. 



3. Results 



3.1. Molecular emission 



The '^CO (1-0), "CO (1-0), and C'^0 (1-0) emissions were observed simultaneously. 
However we did not detect any C'*^0 (1-0) emission in this region above 3cr (0.3 K), which 
indicates low column density in this region. '^CO (1-0) emission was detected in almost the whole 
region even at the position of the Wolf-Rayet star HD 211853. '^CO (1-0) emission was only 
detected in the ring-like region, but not within it. The detailed analysis of the molecular emission 
is described below. 



3.1.1. The systemic velocity of the molecular clumps 

In Figure 2, we present the 1' averaged spectrum of '^CO (1-0) at the position of the 
Wolf-Rayet star. ^^CO (1-0) emission is significantly detected. From the gaussian fit, the peak 
antenna temperature is 0.38+0.02 K, about 8cr level (Icr equal to 0.05 K estimated from the 
baseline). It has a line width of 4.1+0.2 km s~^ and a peak velocity of -43.4+0.01 km s~^ 

The left panel of Figure 3 presents the 14' x 14' averaged spectra of the HI 21 cm line, line 
'^CO (1-0) and line ^^CO (1-0). The right panel of Figure 3 shows the zoomed in averaged spectra 
of the ^^CO (1-0) and ^^CO (1-0) lines. The re are two compone nts in the ^^CO (1-0) and "CO 



(|2010l) . The narrow one peaks at -50 



(1-0) emission, which was also noticed by .Vasquez et al 

km s \ while the other broader one ranges from -48 km s~' to -39 km s~^ 

Are both components related to this region? We argue that the -50 km s~' maybe a foreground 
cold cloud, rather than part of this region based on the following two aspects. 

(1). The integrated emission of ^^CO (1-0) in the first two panels of Figure 5 is concentrated 
in an isolated and compact core (denoted as core 'T") north-west of HD 211853. We did not 
detect this component in the other positions. Additionally, from the P-V diagrams in Figure 6, one 
can se e that the -50 km s~^ component is also not connected to the other component in velocity 



space. 



Vasquez et al 



(|2010|) argued that the ^^CO (1-0) emission between -53.4 and -46.0 km s"^ 
is also associated with the ring nebula. However, their conclusion may be misleading because 
their integrated map is contaminated by the emission between -48 to -46 km s"'. As the channel 
map shows, the emission in their clouds E and F is dominated by the emission between -48 and 
-39 km s ^ rather than the -50 km s"^ component. 

(2). The left panel of Figure 4 presents the spectra of the HI 21 cm line, ^^CO (1-0) and "CO 
(1-0) lines at the peak of core 'T". The right panel of Figure 4 shows the zoomed in spectra of the 
'^CO (1-0) and "CO (1-0) lines. One can see that the HI emission at the peak of core 'T" has a 
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very narrow absorption dip. Together with the narrow line width of the molecular emission, an HI 
Narrow Self- Absorption (HINSA) systern can b e identified, suggesting that core "I" is more likely 



a foreground cold cloud (Li & GoldsmithI 



|2003[) . As shown in Table 2, the kinematic temperature 
of core "I" (14.6+0.3 K) is much lower than that of core "F" (21.7+1.0 K), indicating that core "I" 
is more quiescent than core "F". 

In conclusion, the -50 km s"' molecular component is more likely from a foreground cold 
cloud, rather than a part of the molecular environment surrounding HD 21 1853. However because 
the -50 km s"^ molecular component is very near the ring nebular in the spatial and velocity 
space, it may also interact with the expanding molecular ring, which needs further observations to 
confirm. 

From the Gaussian fit, the averaged ^^CO (1-0) spectrum peaks at -43.6+0.1 km s"\ very 
similar to that of the '^CO (1-0) emission taken from the central Wolf-Rayet star (-43.4+0.01 
km s"'). Thus, we adopt -43.5 km s~^ as the overall systemic velocity of the molecular emissions 
in this region. 

3.1.2. The properties of molecular cores 

Figure 5 displays the channel maps of ^^CO (1-0) emission in contours overlaid on the 
8.3 /im emission in the MSX A band. The middle velocities of each panel are plotted on the 
upper-right comers. In each panel, the integrated velocity interval is 2 km s~'. The contours are 
from 20% to 90% of the peak intensity in each panel. It can be seen that the ^^CO (1-0) emission 
is very clumpy. The blue shifted gas (-46 km s~' panel) distributes in the western and southern 
regions, while the red shifted gas (-40 km s"^ panel) is mainly in the northern and eastern regions, 
indicating a velocity gradient from north-east to south-west. The gas at -42 and -44 km s"^ forms 
a ring-like structure harboring several well separated cores. The ring-like molecular structure 



perfectly coincides with the 8.3 yum emission in the MSX A band. 

Figure 6 presents the P-V cuts along two orientations denoted by the red dashed lines in 
Figure 7. From the P-V diagrams, it is clear that the -50 km s^ component is separated with the 
-43.5 km s"' component. An amazing aspect of the P-V cuts is that the relative velocity of the 
gas increases with the distance to the WR star from near to far, following a "Hubble law" like 
distribution, as denoted by the dashed lines. 

Figure 7 displays the integrated intensity map of ^^CO (1-0) overlaid on the 8.3 yum emission 
in gray scale, which is integrated from -48 to -39 km s"^ The contours are from 10% (~ 8cr) 
to 90% of the peak emission. Eight cores are identified and denoted from "A" to "H". The 
Wolf-Rayet star HD 211853 is drawn as a "cross". IRAS, MSX and AKARI point sources are 
marked with "boxes", "triangles," and "stars", respectively. However, one should keep in mind 
that many of these point sources may be not real; rather, they may correspond to some variations of 
brightness in the infrared emission along the filamentary structures. More observations are needed 
to constrain their properties. For convenience, we still call them "point sources" in the following 
text. These cores form a ring-like structure surrounding HD 21 1853, and they are associated with 
the 8.3 //m emission in the MSX A band. Most of the infrared point sources are distributed in 
the molecular ring. Core "A" is associated with an IRAS source "IRAS 22172-1-5549", while the 
emission peaks of the other cores seem to be separated from the infrared point sources, which will 
be discussed later in Section 4.3. 

Taking into account the effect of beam smearing, the radii of the cores are calculated as 

I 2 

R=-^^-^ — -, wher e A is the projecte d area of each core and 6b is the beam width. Adopting a 



distance of 3 kpc (|Cappa et al. 



120081) . the radii of the cores are calculated and listed in the fourth 
column of Table 1. The average radius is 1 pc. 

Figure 8 presents the spectra of ^^CO (1-0), ^^CO (1-0), and C^^O (1-0) at the emission 
peaks of the cores. The spectra at core "E" are totally blue shifted, while at core "H" the spectra 
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are red-shifted. The peak velocities of the other cores coincide with or slightly deviate from the 
systemic velocity. The results of the Gaussian fit to these spectra are listed in Table 1 . Core "G" 
has the largest line width and antenna temperatures, while core "H" has the smallest ones. The 
line width of ^''CO (1-0) ranges from 1.3 to 2.6 km s"' with a mean value of 2.1 km s '. 

Assuming '^CO (1-0) emission is optically thick and '^CO (1-0) is optically thin, we derived 
the parameters of the cores including the excited temperatures Tg^, column densities of H2, N//,, 
and optic al depth of '^CO (1- 0) lines under the local thermal equilibrium (LTE) assumption 
following . Garden et al.. a991.) . Typical abundance rat ios [H2]/['^CO] = 104 and [i2CO]/['3co]=60 



were used in the calculations (|Deharveng et al. 



20081). The volume density of H2 is calculated 



as n//,=N//2/2R. Then the core mass can be derived as M=|7r • R^ ■ uej ■ niu^ ■ fig, where m//, is 
the mass of a hydrogen molecule and /ig=1.36 is the mean atomic weight of the gas. The derived 
parameters are listed from Columns 3-7 in Table 2. The N//^ ranges from 1.8 x 10^' to 1.4 x 10^^ 
cm"^. The optical depth of '^CO (1-0) is around 0.3. The core masses range from 21 Mq ("H") to 
1061 Mo ("G") with a mean value of 400 Mq. 



We also applied RADEX dVan der Tak et al 



20071), a one-dimensional non-LTE radiative 



transfer code, which uses the escape probability formulation assuming an isothermal and 
homogeneous medium without large-scale velocity fields, to constrain the physical conditions, 
such as density and kinetic temperature, of the cores. In the first running, we explored a range 
of H2 volume densities and temperatures of [10^,10^] cm~^ and [5,50] K. The molecular column 
densities, which are also input parameters required, are fixed according to the H2 volume densities 
as N(X)=n//2-2R-X(x), where X(x) is the relative abundance ratio of the molecule. For ^^CO 
(1-0), X(x)=10"\ and for ^^CO (1-0), X(x) is adopted as 1.67x10 ^ The model parameters were 
selected to satisfy | Tmod-Tohs 1^ 3cr, where Tmod and Tots are the modeled and observed brightness 
temperature of the molecule transition, respectively. As shown in the left panel of Figure 9, 
the kinematic temperatures and H2 volume densities can be well constrained by modeling both 
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the ^^CO (1-0) and '^CO (1-0) lines. Two geometries, static spherical and expanding spherical 
(similar to large velocity gradient (LVG) approximation), were chosen in this running. It can be 
seen that both geometry assumptions give nearly the same results for '^CO (1-0) lines, while for 
^^CO (1-0) lines, static spherical models give slightly higher volume densities than expanding 
spherical at higher kinematic temperatures. However, in the overlapping region, the deviations 
of the two geometries are small enough to be ignored. We narrowed the parameter space in the 
second running assuming the expanding spherical geometry. The range H2 volume densities and 
temperatures in the second running are [2.5x10^,10"^] cm'^ and [8,30] K, respectively. As shown 
in the right panel of Figure 9, the kinematic temperatures and volume densities of H2 can be well 
constrained in the overlapping region of ^^CO (1-0) and '^CO (1-0) tracks. The best physical 
parameters from RADEX are obtained by averaging the parameters in the overlapping region and 
are displayed from Columns 7-14 in Table 2. 

It is interesting to note that the excitation temperatures of ^^CO (1-0) obtained from RADEX 
are perfectly consistent with those calculated under LTE assumption. The excitation temperatures 
of "CO (1-0) are much lower than those of "CO (1-0). The optical depths of "CO (1-0) obtained 
from RADEX are larger than those obtained from LTE conditions. The kinematic temperatures 
are always larger than the excitation temperatures of "CO (1-0). The volume densities of H2 and 
core masses calculated in RADEX roughly agree with those from LTE calculations. 

3.2. Radio continuum emission, HI emission and PAH emission 



( 20101) obtained 



From the radio emission detected in CGPS at 1420 GHz, .Vasquez et al. 
an rms ne=20 cm"^ and M////=1500 M© in this region. As shown in the upper panels of Figure 
10, nine NVSS radio sources (shown in blue dashed contours) are located in this region, and 
most of them are bordered by the 8.3 yum emission ring shown in red solid contours. The nine 
radio sources are not individual HII regions, but are instead high brightness components of the 
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large HII region. They reflect the fluctuation and non- uniformity of the ionized gas distribution 
within the large HII region. The Wolf-Rayet star is associated with NVSS source "1", and the 



Panagia & WalmslevI (119781) . 



molecular core "A" is associated with NVSS source "5". Following 
the electron densities, emission measures and masses of ionized hydrogen of these NVSS sources 
are derived by assuming an electronic temperature of 8000 K. However, those parameters listed 
in Table 3 can only be treated as lower limits due to the missing flux of VLA. NVSS "4", "5" and 
"2" in the yellow dashed box show very high electronic densities and form an ionization front 
interacting with the molecular ring. NVSS "1" and "3" seem to interact with the northern part of 
the molecular ring, mainly with core "G". 

We integrated HI emission from -48 km s"^ to -39 km s"^ and displayed the integrated 
intensity map in gray scale in the upper panel of Figure 10. The HI emission is very diffuse with 
an integrated intensity gradient in the N-S orientation. The HI emission also has a shell-like 
appearance in the southern part, which is associated with the 8.3 /im emission. An HI cavity is 
found within the yellow dashed box, where the ionization front exists. Assuming that the HI 21 
cm line is optically thin and spin temperature Ts » hv/k, the average column density of HI can be 
calculated as Nhi = 1.82 x 10'^ x J Tydv (cm"^). Based on the averaged spectrum of HI shown in 
the left panel of Figure 3, the average column density of HI in the whole region of shell B is about 
9.4 X 10^° cm~^. The total mass of HI is estim ated to be 1 100 M©, which is comparable to the total 



mass of HII (-1500 Mq) dVasquez et al. 



2010 ) 



The ring-like f eature of the 8 . 3 //rn emission in MSX A band was been revealed by 



Cappa et al. 



(120081) : 



Vasquez et al. 



(120101) . It clearly shown as red contours in the upper-panel of 



Figure 10 in this p aper. The 8.3 //m ern ission is dominated by PAH emission and indicates the 



presence of PDRs ([Vasquez et al. 



20101) . IRAC 8 yum emission shown in color scale in the lower 



panel of Fig. 10 also reveals a shell-like structure of PAH emission. 



In Figure 1 1, we analyze the distribution of the normalized intensities of CO (1-0) emission 
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HI emission, 1 .4 GHz radio emission observed by CGPS and PAH emission revealed by MSX A 
band along four orientations centered on the Wolf-Rayet star HD 21 1853. At first glance, one 
can see that ^^CO (1-0) emission and PAH emission encircle a hollow with an angular radius of 
~ 5' in each panel, where their emissions are suppressed but 1.4 GHz radio emission is enhanced, 
indicating that the central region is filled with UV radiation that has photodissociated and ionized 
the molecular and atomic gas. The central Wolf-Rayet star is surrounded by the ionized gas 
and shows strong continuum emission at 8.3 yum. The HI emission always shows distribution 
opposite to the 1.4 GHz radio emission. In other words, the peaks and troughs of HI emission 
always correspond to the troughs and peaks of the 1.4 GHz radio emission, respectively. This 
anti-correlation is more significant in the PA=0° plot. However, we also noticed that the HI 
emission is also high in some regions filled with ionized gas, such as in NVSS 1 and 3. This 
may be caused by the non-uniform distribution of HI gas, or by the projection of the atomic gas 
located at the outer layer of the expanding bubble. The peaks of ^^CO (1-0) emission generally 
are located outside of PAH emission. The fact that PAH emission is located at the interface 
between the ionized and molecular material sug gests that molecular g as is being photodissociated 



by the UV photons, i.e., the presence of PDRs dVasquez et al 



2010 ). At the radius of ~ 5' in 



the PA=135° panel, ^^CO (1-0), HI, 1.4 GHz radio and PAH emissions nearly all reach their 
maximum emissions. This position corresponds to molecula r core "A", vyhich i s a bright Rimmed 



cloud and is forming the youngest star cluster in this region (|Saurin et al. 



2010) 



3.3. Spectral Energy Distributions (SEDs) of YSOs 



Oiu et al. 



(|2008l) identified tens of Young Stellar Objects (YSOs) in the south-east of this 
region. We selected 64 of these sources for analysis. The others are beyond 5' from the Wolf-Rayet 
star HD 211853 and are mostly classified as class II objects. The SEDs of thes e YSOs were the n 



modeled using an online two-dimensional radiative transfer tool developed by 



Robitaille et al 
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(120061 12007|) . As input parameters, the distances of these YSOs and the visual extinction Ay were 
explored in a parameter space of [2.75,3-25] kpc and [2,2.5] mag, respective ly. The adopted inpu t 



A,, is similar to that of HD 21 1853 (2.28 mag) and BD+55°2722 (2.26 mag) dVasquez et al 
The SED models that satisfy the criterion of ;\f^ 



20101) 



^lest < 3 X Hdata, whcrc Hdata is the numbcr of 



data points, were accepted and analyzed. Figure 12 shows the SEDs of ten stars. 



Following 



Grave & Kuman (|2009|) . a weighted mean and standard deviation are derived for 



all the physical parameters of each source from the selected models, with the weights being the 
inverse of the x^ of each model. The fitting results are listed in Table 4. The various columns 
in Table 4 are as follows: Column 1, the name of the star; Columns 2 and 3, the coordinates 
of the star; Column 4, the number of models averaged; Column 5, the weighted ^^; Column 
6, the distance of the star; Column 7, the interstellar extinction; Column 8, the age of the 
star; Column 9, the mass of the star; Column 10, the total luminosity of the star; Column 11, 
the envelope mass; Column 12, the envelope accretion rate; Column 13, the disk mass; and 
Column 14, the disk accretion rate. The evolutionary stages of these YSOs are defined as follows 
(Robitaille et al. 2006): Stage O/I objects are those with M^nvIM^ > \Q~^yr~\ Stage II objects 
are those with MgnvIM^ < lO'^yr'^ and M,usklM^ > 10"^, and Stage III objects are those with 
Menv/M^: < lO'^yr'^ and Mdhk/M^: < 10"^. The identifications of the evolutionary stages are listed 
in the last column of Table 4. "Y14" and "Y20" are badly fitted with;^^^ much larger than 1000, 



indicating that they may not connect to this region. We exclude ther n in the fq l 



Oiu et al. 



owing analysis. 



(l2008h . is more 



The SED of "Y27", which was once classified as a class II YSO by 
likely from the atmosphere of an evolved star rather than from YSOs. It can be well modeled 
by a stellar atmosphere with temperature of 7500 K and log(g) of 0.5. The last panel of Figure 
12 presents its SED. Among the other YSOs, 40% of them are younger than 10^ yr and 44 have 
masses smaller than 2 M©. The youngest and most massive YSOs are found associated with 
molecular core "A". 
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One should be careful with these results for individual young stars due to the limit of bands 
used in fitting. More data, especially at long wavelengths, are required to better constrain the 
properties of these young stars. However, the present results can provide us clues as to the age and 
the mass distribution of these young stars spreading in this region. To investigate the reliability 
of the age and mass distribution of these YSOs in the SED fitting, we fit the SEDs again by 
loosing the parameter space of A,,, which is explored in [0,2.5] mag. Figure 13 compares the 
ages and masses of the YSOs in the two runnings. One can find that they coincide with each 
other very well. We therefore believe that even if the parameters of the individual stars need to be 
further constrained, the trend of the age and mass spatial distribution of these YSOs is reliable in 
statistics. 



3.4. The properties of core A / IRAS 22172+5549 



Molecula r core "A" is ass ociated with IRAS 22172+5549, which is an externally heated 



rimmed cloud (|Oiu et al. 



120081) . The lower panel of Figure 10 displays the 850 /im emission in red 
solid contours from SCUBA, which shows a compact dust core. We collected data from IRAS, 
AKARI and SCUBA, and modeled the SED with wavelength longer than 60 /^m, whose emission 
is mainly from a cool dust envelope, with a simple isothermal gray-body dust model. In the 
optically thin case, the gray-body dust model can be written as (Zhu et al. 2010, and references 
therein): 

where Sy is the continuum emission flux at the frequency v, Mtot is the total mass of the gas and 
dust, Ky is the dust opacity per unit dust mass, g=100 is the density ratio of the gas to dust, D is 
the distance, and By(Trf) is the Planck function with a dust temperature of T^. 

For a typical dust coagulation timescale Tcoag = 10^ yr, in the molecular core "A" 
"ff X Tcoag ~ 10^ cm"^yr< 10^ cm"^yr, indicating that coagulation does not have a considerable 
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effect on the dust opacity in core "A" (|Ossenkopf & Henninglll994|) . Thus, the dust opacity of 
the initial d ust distribution (without dust coagulation) with thin ice mantles can be used in the 



SED fitting dOssenkopf & Henning 



1994|) . The dust opacity at wavelen gth A can be derived as 



K, = .797 X ( 



m-^k 



looo/im 



) ^^^, which is from fitting the data in Table 1 of 



Ossenkopf & Henning 



(119941) . Figure 14 shows the SED of IRAS 22172+5549 overlaid by the best fit curve. The dust 
temperature from the best fit is 26.0+0.6 K, which is larger than the kinematic temperature derived 
from the CO emission, suggesting the dust emission traces a more innermost part than CO gas. 
The total mass of gas and dust traced by dust emission is 102+13 Mq, which is much smaller 
than that traced by CO gas. The bolometric luminosity integrated from 20 to 3000 fxm is 2200 
Lq. The dust temperature and the bolometric luminosity obtained here are consistent with the 
dust temperature (28+5 K) and IR luminosity (2 700 Lq) calculated using only the measured flux 



densities at 60 and 100 //m (|Vasquez et al 



2010 ) 



2010 ). 



An embedded infrared cluster "SBBl" is found in this molecular core (|Saurin et al.l 
From SED modeling, there are found to be associated with this molecular core nine YSOs younger 
than 5 x 10^ yr, three of which are even younger than 5 x 10"* yr, indicating that "SBBl" is a very 
young star cluster. As shown in Figure 12, all nine YSOs have large infrared excesses. The three 
most massive ones are "Y54", "Y55" and "Y56", which have stellar mass > 6 M©. The other 
six sources have stellar mass < 3 Mq. The infall rate of these nine YSOs is found larger than 
~ 1 X 10"^ Mq ■ yr"^ . "Y57" has the largest infall rate of ~ 1 x 10"^ M^ ■ yr'\ All rates indicate 
that young stars are forming in core "A". 
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4. Discussion 

4.1. The geometry of the ring nebular 

The geometry of the ring nebular is important to understanding the interaction between the 
central Wolf-Rayet star and its environment. In this section, we argue that the molecular ring 
surrounding HD 21 1853 is more likely flattened rather than spherical. First, there is less molecular 
emission within the interior of the molecular ring, which excludes the spherical shell hypothesis. 
Second, from the first moment map of ^^CO (1-0) in the left panel of Figure 15, one can find that 
the velocity distribution of the molecular gas is not spherical. The gas is redshifted and blueshifted 
in the north-east and south-west, respectively. But the gas along the NW-SE direction is much less 
shifted, indicating that the molecular gas is distributed in a ring with a thickness much smaller 
than the size of the bubble. The plane in which the molecular ring lies has an inclination angle 
with respect to the sky plane, leading to the velocity gradient of the molecular gas along NE-SW 
direction. Finally, as shown in Figure 2, the spectrum of '^CO (1-0) at the position of the WR 
star only has one velocity component, which excludes the two expanding layer hypothesis. The 



geometry a nd the formation of mo 



in detail by 



Beaumont & Williams 



ecular rings surrounding the expanding bubbles were discussed 



(|2010|) . Li their interpretation, the morphology of molecular 
rings naturally arises if the host molecular clouds are oblate, even sheetlike, with thicknesses of 
a few parsecs. The central massive stars can clear out a cavity in the flattened molecular cloud 
and create the molecular ring surrounding them. This is only in the case of the molecular ring 
surrounding HD 211853. 

A simple geometric model is sketched in the right panel of Figure 15. The axes x and y lie on 
the plane of the sky. The WR star HD 21 1853 is located at the origin. The ring nebular is shown 
as a red circle. The projection of the ring nebular in the sky plane is shown as a blue ellipse. The 
lengths of the semimajor and semiminor axis of the ellipse are denoted as "a" and "b". The angle 
9 is the inclination angle between the plane of the ring nebular and the sky plane. The angle is 
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the position angle (P.A.) measured from the major axis. From this simple geometric model, the 
inclination angle 6 can be inferred from: 

cos0 = - (2) 

a 

The observed velocity, Vobs of a point on the ellipse is given by: 



1 
1 + cos^ 6 ■ cot^ ^ 



Vobs = Vexp -d-sineJ — -— + Vsys (3) 



where Vgyg is the systemic velocity, V^j^p is the expanding velocity, 6=1 when 0° < cp < 180° and 
6=-l when 180° < (f> < 360°. The two dashed ellipses in the left panel of Figure 15 approximately 
mark the boundary of the interior and exterior of the molecular ring. The inner and outer ellipses 
are of 4'. 5 x 3' and 7' x 5' in diameter. Thus the inclination angle between the plane of the ring 
nebular and the sky plane turns out to be ~ 44° - 48°. From the left panel of Figure 15, the 
maximum value of | Vohs - Vsys I is found to be 3.5 km s"^ along the direction of the minor axis, 
indicating an expanding velocity of ~ 5 km s ^ 

As shown in Figure 1, the ring nebular surrounding HD 21 1853 is part (shell B) of the large 
HII region Sh2-132. Although shell B is mainly ionized and created by HD 211853, its interaction 
with the other part (shell A) cannot be ignored. The formation of shell B should be affected by 
the expansion of shell A, especially in the west, where the molecular ring is not closed. In the 
above model, we did not take into account the interaction with shell A, which needs more detailed 
modeling. However, the simple model can roughly depict the ring structure surrounding HD 
21 1853 and explain the spatial and velocity distribution of the molecular gas in the ring nebular. 

4.2. Gravitational stability of the molecular cores 

The gravitational stability of the molecular cores can be investigated by comparing their core 
masses with virial masses and Jeans masses. Assuming that the cloud core is a gravitationally 
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bound isothermal sphere with uniform dens ity and is supported sole ly by random motions, the 



virial mass M,,,> can be calculated following 



Ungerechts et al 



(I2OOOI) as following: 



R AV 



= 2.10x 101— )(- 



■) 



(4) 



Mq pc km s 

where R is the radius of the core and AV is the line width of "CO (1-0). The virial masses are 
listed in the 15th column of Table 2. 

Many factors, including thermal pressure, turbulence, and magnetic field, support the gas 
against gravity collapse in molecula r cores. Taking into account th e thermal and turbulent support. 



the Jeans mass can be expressed as (|Hennebelle & Chabrier 



20081): 



M, 



aj- 






(5) 



where aj is a dimensionless parameter of order unity which takes into account the geometrical 
factor, p is the mass density and Cs,eff is an effective sound speed including turbulent support, 

where C^ is the thermal sound speed and ctnt is the non-thermal one dimensional velocity 
dispersion. The thermal sound speed is related with the kinematic temperature as following: 
Cs = {kTkliXgniH^y^^, where //g=1.36 is the mean atomic weight of the gas. The non-thermal one 
dimensional velocity dispersion <jnt can be calculated as follows: 



O'NT 



crt 



CO 



mnco 



and 



AV 



cruco = 



ylUMT) 



(7) 



(8) 



with mnco being the mass of CO. By introducing an effective kinematic temperature 
Tfff = "" j^^ "^ , equation (5 ) can be rewritten in a form similar to Equation (18) in 



Hennebelle & Chabriej (l2008h . 



M, . ,o.x^,-(^r'«( 



in4 -^ r'^'^o 

10^ cm ^ 



(9) 



20 



where n is the volume density of H2 calculated with RADEX, and ju is the mean molecular weight 
of the gas. The Jeans masses calculated are listed in the 16th column of Table 2. 

The core masses of cores "C", "E", and "H" are found to be much smaller than their viral 
masses and Jeans masses, while the core masses of the other cores are comparable with their 
virial masses and Jeans masses. Keeping this in mind, we did not consider the external pressure 
from their ionized boundary layers in calculating the virial masses and Jeans masses above. 
It seems that these cores are gravitationally stable against collapse without extemal pressure. 
However, the rising extemal pressure surrounding t hem can significantly change their stabilities 



an d induce collapse in them (IThompson et al 



is ([Morgan et al. 



2004): 



2004|) . The external pressure from the ionized layer 



= 2piC^ = AnJe 



(10) 



The effective electron temperature Tg is assumed to be 8000 K. The electron densities n^ in the 
ionized boundary layers surrounding cores "A", "B", and "G" are taken as the val ues of NVSS 



sources 5, 4, and 1, respectively. For the other cores, we assume Ug to be 20 cm ^ (^Vasquez et al 



2010!) . The calculated extemal pressure from the ionized layers can be found in the 18th column 
in Table 2. Molecular pressure inside the clouds can be expressed as: 

"mol 



nT, 



eff 



(11) 



The inferred P,„oi can be found in the last column of Table 2. We can see for all the cores that 
extemal pressure from the surrounding ionized gas is much greater than the molecular pressure, 
indicating that the photoionization-in duced shocks can pe netrate the interiors of the molecular 



cores and compress materials inside ([Morgan et al 



20041) . To explore the stabilities of these 



cores subject to non-ne gligible extemal pressure, we calculated the "pressurized virial mass" as 



([Thompson et al. 



2004) 



9 (AV)^ 
Mp, ^ 5.8 X 10-2-^ — ^ 



(12) 
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where G is the gravitational constant and AV is the line width of ^^CO (1-0) lines. The calculated 
"pressurised virial masses" are listed in the 17th column of Table 2. The molecular cores with 
core masses exceeding their "pressurized virial masses" are unstable against gravitational collapse 



(IThompson et al. 



20041) . Core "C", "E" and "H" have core masses smaller than their "pressurized 
virial masses", indicating they are stable even under external pressure. For the other cores, their 
core masses are much larger than their "pressurized virial masses", suggesting that they are likely 
to be unstable against collapse. It can be seen the high external pressure has a great destabilizing 
effect upon the gravitational stabilities of the cores. 

4.3. The effect of the Wolf -Ray et star on the molecular distribution 

It is expected that the current WR sta rs and their pr e vious main sequence phases should 



greatly influence their associated clouds. iVasquez et al. 



(|2010l) has demonstrated that the 



mechanical energy released by the Wolf-Rayet star HD 21 1853 can shape the ring-like structure 
of the molecular and PAH distributions. 

From Equation (3), one can see that the gases distributed along a direction of a given RA. 
should have the same observed velocity. However, this is not true in the first moment map in the 
left panel of Figure 15. Instead, one can notice that the expanding velocity increases with the 
distance from the origin, especially for the gases along the NE-SW direction. This situation can 
also be found in the P-V diagrams in Figure 6, where the velocity distribution shows a "Hubble 
aw". Such a ve 



( Dentetal. 



ocity gradient was also noticed in the molecular clouds surrounding O stars 
20091). It seems that the expansion of the ionized gas can constantly accelerate the 
molecular gas. However, the acceleration is not significant toward the dense cores as shown in 
Figure 8. The dense molecular cores seem to have slight velocity deviation from the systemic 
velocity, which can be explained by the fact that the expansion of the ionized gas is decelerated by 
the resistance of the dense cores due to their large column densities. HD 21 1853 has a significant 
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effect on its surrounding molecular gas by sweeping, photodissociating, reshaping and possibly 
collecting. 

We also noticed that the infrared point sources in the molecular ring are always separated 
from the molecular cores. This effect is more significant toward cores "C" and "F", where the 
infrared sources are located in front of the molecular core and face the central Wolf-Rayet system, 
indicating that the molecular cores have different expanding velocities with the infrared point 
sources. This separation indicates that stellar wind has different effects on the dust and gas. 
However, such separation maybe caused by the separation of the formed stars from their parent 
clouds during evolution or by inducing new cores by those formed stars. Detailed analysis and 
observations are needed to settle this problem. 

4.4. Sequential star formation in the vicinity of core "A" 

4.4.1. Large-Scale sequential star formation 

YSOs with various ages are drawn with different markers in Figure 16 overlaid on the column 
density of H2 in gray scale and excitation temperature of ^^CO (1-0) in contours. The interesting 
thing is that the age gradient of the YSOs appears in a large-scale. The oldest YSOs with ages 
larger than 2 x 10^ yr are scatted in a large area from the Wolf-Rayet system to the molecular 
ring. The younger ones with ages greater than 5x10^ yr but smaller than 10^ yr are located in 



the midway, while the youngest ones are found associate d with mo 



gradient can also be identified in panel (i) of Figure 2 in 



ecular core "A". Such an age 



Oiu et al 



(l2008h . In their paper, the 



evolutionary stages of the YSOs in this region are identified in color-color diagrams. 

We averaged the ages of the YSOs in each 0'.5 radius bin and present the radial age 
distribution in Figure 17. The projection effect was not taken into account. Since the YSOs are 
mainly distributed in the vicinity of core "A", the projection effect should not greatly affect the 
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trend of the age distribution. We find that the YSOs within 2' from the Wolf-Rayet system are as 
old as 2x10^ yr. Those between 2' and 4' are aged ~ 1.2 x 10^ yr. The youngest stars are beyond 
4' and have an average age of ~ 2 x 10^ yr. It is clear that the nearer YSOs are older than the 
farther ones, indicating a stellar age gradient. The stellar age distribution suggests that at least 
three generations of stars have been forming in this region. The oldest generation is closer to 
the Wolf-Rayet system, while the youngest one is l ocated in the mole cular ring. In addition, the 



three star clusters in this region have different ages (jSaurin et al 



2010 ). The "Teutsch 127", which 



includes the 0-type star BD+55°2722 and near the Wolf-Rayet star HD 21 1853, is the oldest one 
(~5 Myr). The "SBB 2" is as old as ~2 Myr. The age of the youngest cluster "SEE 1" is about 1 
Myr. The age gradient of these three star cluster also indicates sequential star formation on a large 
scale. Such "relay star formation" hints at trigged star for mation due to the expansion of bubbles 



created by stellar wind of massive stars (jChen et al 



2007) 



4.4.2. Small-Scale sequential star formation 



As shown in Figure 16, a small-scale age gradient is found to be associated with molecular 
core "A," a bright rimmed cloud. The older and less massive YSOs embrace the molecular core, 
while the youngest and more massive stars reside deeper in the core. This picture is in ac cordanc e 



with th e so-called "radiation-driven implosion" (RDI) process as depicted in Fig ure 1 of 



Qgura 



(120061). which is propose d for triggered star formation in a bright rimmed cloud (Eertoldi 



1989 



Eertoldi & McKee 



19901) . In the model, the surface layer of the cloud is ionized by the UV 



photons from massive stars and then the cloud is compressed and collapsed due to the shock 
generated from the ionization front. As discussed in Section 4.2, core "A" is likely to collapse due 
to extemal pressure from the ionized boundary layer. Recent simulation suggests that there is a 
range of ionizing fl uxes in the RDI model, that trigger star formation in the bright rimmed cloud 



(lEisbas et al. 



201 1|) . This range is about lO^cm ^s ' < ^lyc ^ 3 x 10^'cm ^s ^ If the ionizing 
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flux is larger than lO^cm ^s \ the shock front preceding the D -type ionization fro nt can propagate 



into the cloud, compress it, and trigger star formation therein (|Bisbas et al 



ionizing flux can rapidly disperse the cloud without triggering star formation (IBisbas et al. 



The ionizing fluxes surrounding core "A" can be estimated as ([Morgan et al. 

^LyC 



20 Uh. However large 



2004 ) 



2011 



-0.35,. 0.1/1-2 



i.24xio"'5vr;''y^'-'r 



(13) 



where 5y is the integrated radio flux in mJy, T^ is the effective electron temperature of the ionized 
gas in K, v is the frequency in GHz and 9 is the angular diameter over which the emission is 
integrated in arcseconds. Assuming T^ = 8000 K, the ionizing flux in NVSS 5 is inferred as 
1.1 X 10'°cm"^s \ which is taken as the value in the ionizing boundary layer surrounding core 
"A". This value i s moderate, indicati ng that triggered star formation can take place in core "A". In 



the simulation of 



Bisbas et al 



(120111) . however, the RDI model favors the formation of low mass 
stars in the clouds, which contradicts the situation in core "A", where YSOs with masses larger 
than 6 Mq are found. We noticed that the initial cloud mass in their simulation is much smaller 
than the mass of the core "A", which may be responsible for the contradiction. Additionally, other 
processes such as "collect and collapse" may also play important roles in star formation in core 
"A", which will be discussed in the next section. 



4.5. Collect and Collapse scenario 



The ring-like structure of dust and gas distribution and the existence of PDR in- 



dicate that the star formation in this region favors t 



(IDeharveng. Zavagno. & Caplan 



2005 



Vasquez et al. 



le "co llect and collapse" scenario 



2010|) . The molecular ring is composed of 



well separated, dense and massive cores, most of which are collapsing. Such regular spacing 
rule of the cores could not be pre-existing clumps but are formed by extemal triggering, which is 



a strong evid ence for the "collect and collapse" process (IDeharveng. Zavagno. & Caplan 



2005 



Ogura 



2010|) . This situation contradicts the RDI model, in which the dense clumps into which 
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the radiatively induced shocks drive are pre-existing. Additionally, the large scale sequential star 
formation disc ussed in Section 4.4.1 also favors the " collect and collapse" process rather than the 



RDI" model (lElmegreen & Ladt 



1977 



Ogural 



2010h 



However, as discussed in Section 4.4.2, the bright rimmed nebular core "A" also seems to 
favor the "RDI" model. The "collect and collapse" process and the "Radiation Driven Implosion" 
process cannot take place at the same time and place. To resolve this contradiction, we propose 
that these two processes may function at different times in this region. First, neutral materials are 
collected and accumulated into a dense shell between the ionization front (IF) and shock front 
(SF) generated by the expanding HII region. Second, the compressed shocked layer becomes 
gravitationally unstable and fragments into regularly spaced massive clumps on a long timescale. 
These two steps mirror the "collect and collapse" process. Third, the clumps are illuminated and 
eroded by the ionization front, leading to the formation of bright rims and/or cometary globules 
that collapse to form stars due to the external pressure as discussed in Section 4.2. This third step 
is like the "RDI" process. Besides, the ionization front illuminating the bright rimmed nebular 
core "A" is mor e likely generated from young cluster "SBB 2" rather than from the central 



Wolf-Rayet star (|Saurin et al. 



2010h . Relatively speaking, core "A" is pre-existing compared with 



the ionization front generated from the older cluster. In addition, a possible bow- sho ck seems to 
have been generated in the cluster "SBB 2" and seems to have an impact on "SBB 1" (iSaurin et al 
20101). We therefore suggest that the star formation in core "A" may also be affected by the older 
cluster "SBB 2" through a process like "RDI", leading to the formation of the youngest embedded 
star cluster "SBB 1". 



4.6. Comparison with the other star forming regions 

It is necessary to compare the characteristics of the molecular cores in this region with those 
of the other star forming regions. The line width of ^^CO (1-0) ranges from 1.3 to 2.6 km s~' with 



a mean value of 2.1 km s ^ , which i s 



regions (~2 km s ) (ILiu et al. 



2011 
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very similar to that of the intermediate-mass star forming 



Sun et al. 



2006h. It is lar ger than that of the low-mass star 



1983|) . but much smaller than that of the 



forming regions (~1.3 km s ^) (|Myers. Linke & Benson 

high-mass star formin g regions associated with IRAS sou rces (> 3 km s~ ^) (|Wang et al. 



Wu et al. 



2003 



200 Ih or methanal masers (~4.5 km s ) (jLiu et al 



20091 : 



2010h . Judging from the 



line width of ^''CO (1-0) in this region, the molecular cores where the triggered star formations 
possibly take place favor formation of intermediate-mass YSOs, which is confirmed by SED 
modeling of the YSOs in molecular core A. As discussed in Section 4.4.2, core A is suggested to 
suffer an RDI process. A tendency toward HAeBe stars f ormin g deeper in the cloud is expected 



in RDI model dOgura 



20061 : 



Chen et al. 



2007 ) 



Liu et al. 



(|201l|) found that the regions harboring 



HAeBe stars younger than 10^ yr have an average column density of ~ 5 x 10^' cm~^, slightly 
smaller than that in most cores of this region. They also have an average excitation temperature 
of ~16 K in those HAeBe regions younger than 10^ yr, which is slightly lower than that of the 
molecular cores of this region (~20 K), indicating the influence of external heating. 



5. Summary 

We studied the environment surrounding Wolf-Rayet star HD 21 1853 in molecular, infrared, 
radio and HI emission. The SEDs of the YSOs that are located in the south-east part of this region 
are modeled with an online YSO SED fitting tool. The triggered star formation is discussed. The 
main findings of this paper are as follows: 

(1). The molecular emission has two components. One peaks at -50 km s ', which is likely 
from a foreground cold cloud. The other one has a systemic velocity of -43.5 km s ', whose 
'^CO (1-0) emission forms a ring-like structure with several well separated molecular cores. The 
line width of '^CO (1-0) ranges from 1.3 to 2.6 km s"' with a mean value of 2.1 km s^\ which is 
very similar to that of intermediate-mass star forming regions. 
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(2). The properties of the molecular cores are analyzed under the LTE assumption as well 
as with non-LTE models in RADEX. Both LTE and non-LTE analysis give the same excitation 
temperatures for ^^CO (1-0) (~20 K). The excitation temperatures of ^^CO (1-0) calculated with 
RADEX are much smaller than those of ^^CO (1-0). The core masses calculated under LTE and 
non-LTE are similar, ranging from ~ 10^ to ~ 10^ M©. The volume densities of H2 in this cores 
are ~ lO-' cm"^. 

(3). The molecular emission is associated with 8.3/im emission in the MSX A band. We 
also witnessed structures like cavities and shells in HI emission. Radio emission at 1.4 GHz is 
embraced by the molecular and dust ring. Ionization fronts easily can be identified between the 
radio and molecular emission especially in the south-east part. 

(4). We find that the molecular gas surrounding HD 21 1853 is not spherically expanding. 
The molecular gas distributes in a ring-like structure with a thickness much smaller than the size 
of the bubble. The inclination angle between the plane of the ring nebular and the sky plane turns 
out to be ~ 44° - 48°. The expanding velocity of the ring nebular is about ~ 5 km s"^. The 
expanding velocity seems to increase with the distance from the WR star. 

(5). The molecular cores are well separated along the ring. It seems that these cores are 
gravitationally stable against collapse without external pressure. However the high external 
pressure from the ionized gas has a great destabilizing effect upon the gravitational stabilities 
of the cores. Core "C", "E" and "H" have core masses smaller than their "pressurized virial 
masses", indicating that they are stable even under external pressure. For the other cores, their 
core masses are much larger than their "pressurized virial masses", suggesting that they are likely 
to be unstable against collapse. 

(5). From the SED modeling of the YSOs located in the south-east of HD 21 1853, a stellar 
age gradient is found across a large scale in space from the Wolf-Rayet star to the molecular ring. 
The presence of PDR, the fragmentation of the molecular ring, the collapse of the cores, and the 
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large scale sequential star formation indicate that the "collect and collapse" process can work in 
this region. 

(6). Molecular core "A" is forming a star cluster as young as ~ 10^ yr. Core "A" is embraced 
by the ionization front and older less massive YSOs. More massive YSOs are forming deeper 
in core A. The ionizing flux from the ionized boundary layer surrounding core "A" is moderate 
enough to generate a leading shock, which propagates into the cloud, compresses it, and triggers 
star formation therein. This is very similar to the picture depicted by the "Radiation-Driven 
Implosion" models. 

In conclusion, star formation in the molecular ring is induced by the Wolf-Rayet star HD 
21 1853 through the interaction between the expanding ionizing gas and the ISM. The evolution 
of the molecular ring and the triggered star formation may be aff'ected by both the "collect and 
collapse" and the "Radiation-Driven Implosion" processes. Generally speaking, the materials in 
the host molecular cloud surrounding the WR star are collected into a molecular ring. Then the 
ring fragments into well separated cores that favor the "collect and collapse" process. The cores 
are unstable due to the external pressure and collapse inside to form stars. The ionization front 
proceeded by a shock front generated from the older star cluster (like "SBB 2") can propagate 
into the cloud (core "A") and compress it and trigger star formations therein as depicted by the 
"radiation-driven implosion" process. 
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Fig. 1 . — Composite color image in this region: DSS2-R optical emission is blue, 1420 GHz radio 
emission detected in CGPS is green, and 8.3 /um emission detected at MSX A band is red. The 
larger "star" shows the location of the WR star HD 211853 and the small "star" represents the 
position of the massive 0-type star BD+55°2722. The yellow dashed "box" shows the region of 
shell B, on which we focus in this paper. 
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Fig. 2. — Spectrum of ^^CO (1-0) at the position of the Wolf-Rayet star HD 21 1853 averaged over 
r. The green solid line is the Gaussian fit curve. 
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Fig. 3. — Left: Spectraof ^^CO (1-0), ^^CO (1-0) and HI averaged overan area of 14'x 14' centered 
at the Wolf-Rayet star HD 211853. Right: Zoomed in averaged spectra of ^^CO (1-0), and "CO 
(1-0). 
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Fig. 4. — Left: Spectra of ^^CO (1-0), "CO (1-0) and HI taken from the emission peak of the 
component at -50 km/s (core "I"). Right: Zoomed in spectra of '^CO (1-0), and "CO (1-0) at core 
"I". 
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Fig. 5. — Channel maps of '^CO (1-0). The middle velocities are plotted on the upper-right comers 
of each panel. The integrated velocity interval in each panel is 2 km s'K The contours are from 20% 
to 90% of the peak intensity (4.6 K-km s"'). The background image shows the 8.3 fim emission 
detected at the MSX A band. The position of core 'T" is marked in the -50 km s"' panel. 
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Fig. 6. — P-V diagrams of ^^CO (1-0) cut along the orientations denoted in Figure 7. The contours 
are from 0.6 K (3cr) in steps of 1 K. 
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Fig. 7. — ^^CO (1-0) integrated intensity map overlayed on MSX A band emission. The contours 
are from 10% to 90% of the peak (7.8 K- km s"^). "Triangles" represent MSX point sources, 
"boxes" IRAS point sources, and "stars" AKARI point sources. The Wolf-Rayet star is drawn as 
"cross". The molecular cores are denoted from "A" to "H". The horizon and vertical dashed lines 
represent the directions of the P-V cuts in Figure 6. 
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Fig. 8. — Spectra of '^CO (1-0) are red, those of "CO (1-0) are green and those of C^^O (1-0) are 
blue. The core names are labeled on the upper-left comers of each panel. The dashed black lines 
indicate the systemic velocity of -43.5 km s"^. 
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Fig. 9.— Non-LTE model fitting of ^^CO (1-0) and '^CO (1-0) at the emission peak of core "A". 
Left: parameters run through T/.,„ of [5,50] K and n//, of [10^,10^] cm~^ in log space. Right: 
parameters run through Tyt,„ of [8,30] K and Uh^ of [2.5x10^,10"^] cm"^ in log space. 
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Fig. 10. — Upper: NVSS 1.4 GHz emission is shown in green dashed contours (0.45 mJy (1 
cr) X 3,6,9,...), MSX A band emission is shown in red solid contours (from 10% to 90%), and 
HI emission integrated from -48 km s"^ to -39 km s~^ is in gray scale. Lower: NVSS 1.4 GHz 
emission is shown in green dashed contours, SCUBA 850 //m emission is in red solid contours 
(from 10% to 90%) and IRAC 8 fim emission is in color scale. 
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Fig. 11. — Normalized intensity distributions of CO, HI, 1420 MHz radio emission detected in 
CGPS and 8.3 //m emission detected in the MSX A band along four orientations centered on the 
Wolf-Rayet star. 
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Fig. 13. — Comparison of the stellar mass and age of the YSOs in two runnings of the SED 
modeling. Left: Age comparison; Right: Mass comparison. The solid lines represent the linear 
fits. Age_l and Mass_l are the stellar age and mass obtained in the first running by exploring Ay 
in [2,2.5] mag. Age_2 and Mass_2 are the stellar age and mass obtained in the second running by 
exploring Av in [0,2.5] mag. 
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Fig. 14.— SED of molecular core "A" (IRAS 22172+5549). The solid curve represents the best 
fit. 
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Fig. 15. — Left: first moment map of '^CO (1-0). Two ellipses are drawn for reference, to mark 
the boundaries of the interior and exterior of the molecular ring. Right: Coordinate system for 
describing the orientation in the sky of the ring nebular. The axes x and y lie on the plane of 
the sky. The WR star HD 211853 is located at the origin. The ring nebular is shown as a red 
circle. The projection of the ring nebular in the sky plane is shown as a blue ellipse. The lengths 
of the semimajor and semiminor axis of the ellipse are denoted as "a" and "b". The angle Q is 
the inclination angle between the plane of the ring nebular and the sky plane. The angle is the 
position angle measured from the major axis. 
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Fig. 16. — Tex in contours overlaid on the column density of H2 in gray scale. The contour levels 
are from 20% to 90% of the peak (24 K). The young stars in different age intervals are denoted 
as different markers on the plot as shown in the lower- left comer. The upper age of each interval 
is shown in the lower-left comer. The nine YSOs associated with core "A" are shown as color 
markers. The three most massive YSOs are shown as markers with red edges. The other six less 
massive ones are presented as markers with blue edges. 
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Fig. 17. — Distribution of 0'.5 bin averaged ages of YSOs located between HD 211853 and core 
"A" as a function of radii from the WR star. The errorbars represent the standard deviation in each 
bin. 
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Table 1 . Observed parameters of CO cores 



Core 


RA (J2()0()) 


DEC (J2000) 


R" 


v,.„ 


Ta("CO) 


AV("C0) 


Ta('^CO) 


Av('2C0) 




(hm s) 


(O , „) 


(pc) 


(kms-') 


(K) 


(kms-') 


(K) 


(kms-') 


A 


22:19:09.4 


56:04:46.5 


1.0 


-43.8(0.3) 


2.6(0.1) 


2.4(0.2) 


9.2(0.3) 


3.4(0.2) 


B 


22:18:38.4 


56:03:25.4 


1.4 


-44.2(0.1) 


2.2(0.1) 


2.2(0.1) 


7.5(0.1) 


2.7(0.2) 


C 


22:18:20.1 


56:06:00 


0.9 


-43.7(0.2) 


1.2(0.1) 


2.5(0.1) 


8.5(0.2) 


2.6(0.1) 


D 


22:18:03.1 


56:08:07.6 


1.0 


-43.3(0.2) 


2.5(0.1) 


1.7(0.2) 


7.6(0.3) 


2.6(0.2) 


E 


22:18:03.1 


56:05:42.7 


0.8 


-46.1(0.2) 


1.8(0.2) 


2.0(0.2) 


7.9(0.2) 


2.5(0.2) 


F 


22:18:08.8 


56:11:17.3 


1.1 


-42.7(0.1) 


2.1(0.2) 


1.9(0.2) 


7.7(0.3) 


2.8(0.2) 


G 


22:18:39.1 


56:11:17.3 


1.3 


-43.3(0.1) 


2.9(0.1) 


2.6(0.2) 


10.1(0.2) 


3.5(0.1) 


H 


22:19:26.6 


56:11:17.3 


0.5 


-42.0(0.1) 


1.0(0.1) 


1.3(0.2) 


3.7(0.2) 


2.1(0.1) 


I 


22:18:17.6 


56:12:23.4 


1.0 


-50.1(0.1) 


1.4(0.1) 


1.1(0.1) 


4.9(0.1) 


1.6(0.1) 



"Taking the distance as 3 kpc 



Table 2. Derived parameters of CO cores 



LTE non-LTE 



Core r (I'CO) T„ N^^ n«, M™„ T«„ N„, n„, T„ ('^CO) r ('^CO) T„ ("CO) r ("CO) M„„ M„, My M„, P„,/k P„,„,/k 

(K) (lO^i cm-2) (10^ cm-3) (Mq) (K) (10^' cm-^) (10^ cm"') (K) (K) (Mq) (Mq) (Mq) (Mq) (10^ cm^^K) (10^ cm-^K) 

501 24.0(0.8) 9.0(0.4) 1.4(0.1) 21.7(0.5) 10.6(0.8) 15.6(0.5) 

661 20.4(0.3) 7.8(0.3) 0.9 18.5(0.2) 15.6(0.7) 10.8(0.2) 

193 26.5(1.1) 3.9(0.2) 0.7 20.5(0.4) 7.1(0.5) 9.9(0.4) 

322 20.4(0.9) 6.9(0.3) 1.1(0.1) 18.6(0.5) 14.3(1.2) 12.2(0.3) 

127 22.1(0.8) 5.0(0.5) 1.0(0.1) 19.3(0.4) 10.3(1.1) 11.6(0.6) 

356 21.7(1.0) 6.1(0.6) 0.9(0.1) 18.9(0.5) 11.5(1.2) 11.1(0.5) 

1061 26.2(0.6) 11.3(0.4) 1.4(0.1) 23.6(0.4) 11.2(0.6) 16.7(0.4) 0.6 877(36) 710 687 631 1.1 6.0 | 

21 11.7(0.6) 2.4(0.2) 0.8(0.1) 10.8(0.4) 15.8(1.9) 6.5(0.2) 

99 14.6(0.3) 3.2(0.2) 0.5 13.2(0.2) 19.9(1.2) 7.0(0.2) 1.3 149(10) 231 103 27 0.6 0.4 I 



A 


0.3 


22 


10.9 


1.8 


B 


0.3 


18.5 


7.3 


0.8 


C 


0.2 


20.5 


5.2 


0.9 


D 


0.3 


17.8 


7.0 


1.1 


E 


0.2 


19 


4.3 


0.9 


F 


0.3 


18.8 


6.4 


0.9 


G 


0.3 


23.6 


13.7 


1.7 


H 


0.3 


11.2 


1.8 


0.6 


I 


0.3 


13.2 


2.2 


0.3 



0.5 


410(22) 


504 


544 


199 


5.8 


5.1 


0.9 


701(30) 


647 


523 


174 


3.8 


2.8 


0.5 


145(10) 


473 


870 


703 


0.6 


2.8 


0.8 


319(16) 


357 


232 


150 


0.6 


2.1 


0.6 


149(16) 


336 


383 


288 


0.6 


2.6 


0.8 


339(323) 


439 


350 


235 


0.6 


2.1 


0.6 


877(36) 


710 


687 


631 


1.1 


6.0 


0.9 


27(3) 


137 


121 


51 


0.6 


0.9 


1.3 


149(10) 


231 


103 


27 


0.6 


0.4 



o 



51 



Table 3. Parameters of NVSS radio sources 



Number 



(h m s) 



3EC (J20()0) 


Sl.4 


a 


b 


Or 


n. 


EM 


Mh;/ 


r ' ") 


(mJy) 


(") 


(") 


(') 


(cm"3) 


(pc cm-«) 


(Mo) 


+56:08:24.1 


32.9±4.1 


125.9 


110.0 


1.0 


34 


1322 


1.6 


+56:05:34.6 


4.0±0.5 


<80.5 


<43.7 


<0.5 


<33 


<633 


<0.2 


+56:09:54.6 


20.3±2.5 


83.7 


65.7 


0.6 


54 


2055 


0.6 


+56:04:04.0 


103.4±4.1 


103.0 


56.1 


0.6 


117 


9961 


1.5 


+56:04:42.5 


87.6±3.4 


98.6 


29.3 


0.4 


181 


16879 


0.8 


+56:13:07.7 


10.0±1.4 


71.7 


<50.6 


<0.5 


<52 


<1534 


<0.3 


+56:00:45.0 


18.0±1.8 


52.4 


47.5 


0.4 


92 


4026 


0.3 


+56:12:02.9 


15.2±1.0 


27.5 


<31.8 


<0.2 


<185 


<9675 


<0.1 


+56:00:45.0 


38.2±2.7 


83.2 


64.5 


0.6 


75 


3962 


0.9 



1 221840+560824 22:18:40.77 

2 221846+560534 22:18:46.86 

3 221856+560954 22:18:56.55 

4 221852+560404 22:18:52.33 

5 221905+560442 22:19:05.13 

6 221826+561307 22:18:26.97 

7 221900+560045 22:19:00.72 

8 221929+561202 22:19:29.67 

9 221912+560045 22:19:12.79 



Table 4. The SED fitting results of the YSOs . 



Name RA (J2000) DEC (J200()) Number 
(h m s) C ") 



log(d) 


A. 


log(Age) 


M, 


L, 


og(kpc) 


(mag) 


log(yr) 


Mo 


Lo 



log(Me„v) log(M„i.) log(Mrf,jt) log(Mrf,,,t) Stage 
log(Mo) log(Moyr') log(Mo) logCMoyr') 



Yl 22:18:45.18 +56:04:49.3 35 

Y2 22:18:45.73 +56:03:52.6 100 

Y3 22:18:46.31 +56:05:03.4 306 

Y4 22:18:47.26 +56:05:58.4 37 

Y5 22:18:47.54 +56:05:11.2 107 

Y6 22:18:48.01 +56:05:02.9 1 

Y7 22:18:49.37 +56:05:46.5 1 

Y8 22:18:49.40 +56:04:31.0 1172 

Y9 22:18:49.86 +56:06:18.5 879 

YIO 22:18:50.13 +56:03:40.8 24 

Yll 22:18:50.78 +56:05:06.8 835 

Y12 22:18:52.22 +56:06:13.7 306 

Y13 22:18:52.49 +56:06:48.8 550 

Y\r 22:18:52.69 +56:06:05.1 1 

Y15 22:18:53.62 +56:06:24.8 465 

Y16 22:18:54.35 +56:04:59.6 81 

Y17 22:18:55.39 +56:05:12.9 259 

Y18 22:18:55.67 +56:05:16.3 44 

Y19 22:18:55.67 +56:04:55.4 2584 

Y20'' 22:18:56.14 +56:05:48.6 1 

Y21 22:18:56.27 +56:05:32.3 894 

Y22 22:18:57.18 +56:05:02.0 1 

Y23 22:18:57.62 +56:07:07.0 43 



20.8(1.5) 0.48 2.33(0.04) 6.07(0.08) 1.54(0.20) 8(2) -1.36(0.42) -5.72(0.35) 

12.3(1.0) 0.48(0.01) 2.29(0.06) 5.43(0.45) 1.18(0.26) 20(9) -0.62(0.34) -5.02(0.15) 

10.5(0.8) 0.47(0.01) 2.26(0.07) 5.88(0.24) 0.44(0.14) 1 -2.00(0.32) -5.93(0.31) 

12.6(1.7) 0.47(0.01) 2.28(0.06) 5.40(0.10) 2.08(0.36) 37(18) -0.56(0.26) -4.71(0.34) 



8.0(1.1) 0.47(0.01) 2.27(0.08) 5.95(0.33) 0.24(0.07) 1 

58.1 0.46 2.00 5.18 1.98 31 

130.1 0.46 2.46 5.54 1.30 17 

5.6(1.4) 0.47(0.01) 2.28(0.10) 5.86(0.42) 0.47(0.28) 4(7) 

14.3(0.6) 0.48(0.01) 2.25(0.06) 6.26(0.12) 1.04(0.16) 2(1) 



-1.26(0.62) -5.48(0.37) 

-1.08 -4.74 

-1.11 -5.75 

-0.61(1.36) -4.94(1.00) 

-2.05(0.65) -6.14(0.58) 



11.7(1.8) 0.47(0.01) 2.20(0.06) 6.29(0.18) 2.41(0.35) 28(8) -0.73(0.28) -4.99(0.24) 



5.1(1.4) 0.48(0.01) 2.26(0.10) 6.12(0.29) 0.52(0.19) 1(1) 

7.4(1.1) 0.48(0.01) 2.25(0.08) 5.78(0.40) 0.55(0.25) 4(4) 

12.3(1.1) 0.48(0.01) 2.28(0.06) 6.06(0.16) 0.71(0.15) 2(3) 

30154.2 0.51 2.00 6.09 3.39 36 

15.0(1.2) 0.48(0.01) 2.34(0.05) 5.92(0.14) 0.64(0.12) 2 

6.9(1.3) 0.48(0.01) 2.29(0.08) 5.26(0.70) 0.42(0.14) 4(2) 

5.1(1.5) 0.48(0.01) 2.26(0.10) 6.12(0.30) 0.95(0.30) 3(2) 



-1.77(0.60) 
-1.09(0.50) 
-1.34(1.35) 

-2.98 
-1.46(0.58) 



-5.79(0.54) 
-5.24(0.35) 
-5.43(1.83) 

-8.68 
-5.63(0.42) 



-0.70(0.24) -4.96(0.20) 
-1.06(0.78) -5.30(0.57) 



7.6(1.3) 0.47(0.01) 2.22(0.07) 6.08(0.24) 1.53(0.42) 11(7) -0.55(0.44) -4.91(0.39) 



1.7(1.9) 0.47(0.02) 2.22(0.16) 5.98(0.57) 0.36(0.28) 1(1) 

2131.5 0.48 2.00(0.01) 5.61 4.42 63 

7.1(1.2) 0.48(0.01) 2.27(0.08) 6.19(0.23) 0.83(0.25) 2(1) 

162.5 0.48 2.00 6.38 3.13 128 

46.7(0.8) 0.48 2.32(0.03) 5.84(0.08) 1.19(0.10) 6(1) 



-1.89(2.11) -5.87(0.83) 
-1.10 -5.32 

-1.43(0.59) -5.56(0.48) 
-2.63 

-1.22(0.14) -5.34(0.20) 



-1.90(0.11) 
-1.76(0.22) 
-2.83(0.29) 
-1.95(0.27) 
-2.52(0.20) 

-1.35 

-1.20 
-2.36(0.54) 
-2.59(0.31) 
-1.49(0.17) 
-2.66(0.33) 
-2.42(0.24) 
-2.37(0.28) 

-2.20 
-2.33(0.20) 
-2.14(0.29) 
-2.40(0.30) 
-2.07(0.19) 
-2.91(0.89) 

-2.86 
-2.46(0.31) 

-7.27 
-1.99(0.10) 



-7.13(0.11) 
-6.20(0.41) 
-8.56(0.39) 
-6.68(0.31) 
-7.64(0.42) 

-6.75 

-5.88 
-7.25(0.97) 
-8.22(0.82) 
-6.69(0.31) 
-8.01(0.64) 
-7.64(0.25) 
-7.86(1.87) 

-7.59 
-8.14(0.24) 
-6.00(0.47) 
-7.81(0.36) 
-7.37(0.30) 
-8.45(4.77) 

-9.64 
-8.06(0.39) 

-13.27 
-7.54(0.09) 



O/I 
O/I 
O/I 
O/I 
O/I 
O/I 
O/I 
O/I 
III 
O/I 
O/I 
O/I 
O/I 

111 

O/I 
O/I 
O/I 
O/I 
O/I 
O/I 
O/I 

111 

O/I 



I 



Table 4 — Continued 



Name 


RA (J2000) 


DEC (J2000) 


Number 


7 

X 


log(d) 


Av 


log(Age) 


M, 


L. 


log(M„„,) 


log(M,„,) 


log(Md,jt) 


log(Mrf,jt) 


Stage 






(hms) 


(° ' ") 






log(kpc) 


(mag) 


log(yr) 


Mo 


Lo 


log(Mo) 


log(Moyr-l) 


log(Mo) 


log(Moyr-') 






Y24 


22:18:57.71 


+56:04:27.5 


1084 


4.5(1.7) 


0.48(0.01) 


2.28(0.10) 


6.11(0.35) 


0.45(0.24) 


2(5) 


-1.05(2.14) 


-5.28(1.84) 


-2.50(0.36) 


-7.52(1.26) 


o/i 




Y25 


22:18:57.82 


+56:05:52.6 


270 


16.0(0.8) 


0.47(0.01) 


2.27(0.06) 


5.34(0.47) 


0.44(0.18) 


5(5) 


-0.42(0.65) 


-4.90(0.50) 


-2.29(0.39) 


-6.84(1.40) 


o/i 




Y26 


22:18:58.38 


+56:04:23.3 


24 


15.1(1.4) 


0.48 


2.34(0.06) 


6.88(0.04) 


2.24(0.28) 


46(29) 


0.16(0.39) 


-6.38(0.39) 


-2.21(0.25) 


-7.46(0.37) 


III 




Y27b 


22:18:58.68 


+56:07:23.6 


1 


19476.0 


0.51 


2.50 


6.11 


14.60 


19200 


-7.86 




-8.04 


-9.44 


III 




¥28 


22:18:58.71 


+56:04:54.5 


53 


5.3(1.6) 


0.48(0.01) 


2.24(0.09) 


5.50(0.76) 


0.55(0.39) 


6(8) 


-0.55(0.95) 


-4.99(0.75) 


-2.41(0.67) 


-7.32(0.35) 


O/I 




Y29 


22:18:58.76 


+56:05:18.8 


1106 


4.6(1.6) 


0.47(0.01) 


2.26(0.10) 


5.55(0.38) 


0.26(0.15) 


1(1) 


-1.61(0.67) 


-5.60(0.62) 


-2.85(0.42) 


-8.17(1.60) 


O/I 




Y30 


22:18:59.21 


+56:05:00.6 


62 


2.9(2.1) 


0.47(0.01) 


2.27(0.10) 


5.46(0.77) 


0.56(0.25) 


3(4) 


-1.07(0.43) 


-5.23(0.51) 


-2.19(0.34) 


-7.95(0.56) 


O/I 




¥31 


22:18:59.72 


+56:05:13.9 


1 


128.7 


0.49 


2.02 


4.58 


1.67 


33 


-0.33 


-4.82 


-1.24 


-6.01 


O/I 


1 


Y32 


22:19:00.08 


+56:03:14.8 


314 


18.4(0.7) 


0.48(0.01) 


2.31(0.04) 


6.88(0.02) 


2.36(0.06) 


31(3) 


-3.63(0.21) 




-3.06(0.29) 


-7.99(0.52) 


III 


1 


Y33 


22:19:00.18 


+56:05:34.7 


55 


6.0(1.5) 


0.47(0.01) 


2.22(0.08) 


6.01(0.36) 


0.32(0.11) 


2(1) 


-1.17(0.39) 


-5.29(0.40) 


-2.35(0.21) 


-7.26(0.43) 


O/I 


1 


¥34 


22:19:00.26 


+56:06:41.4 


178 


4.4(1.6) 


0.48(0.01) 


2.22(0.09) 


5.59(0.57) 


0.71(0.35) 


6(8) 


-0.72(0.60) 


-5.03(0.42) 


-2.36(0.35) 


-7.50(0.56) 


O/I 




Y35 


22:19:02.01 


+56:07:08.5 


7 


22.0(1.2) 


0.47 


2.15(0.04) 


6.43(0.09) 


2.30(0.11) 


12(2) 


-1.65(0.23) 


-5.96(0.24) 


-1.85(0.19) 


-7.75(0.24) 


O/I 




Y36 


22:19:02.01 


+56:08:00.7 


1391 


7.8(1.1) 


0.48(0.01) 


2.27(0.08) 


5.85(0.33) 


0.34(0.16) 


1(2) 


-1.32(1.17) 


-5.48(1.20) 


-2.64(0.34) 


-7.68(1.05) 


O/I 




¥37 


22:19:02.18 


+56:07:25.8 


1174 


3.7(1.7) 


0.48(0.01) 


2.28(0.11) 


5.98(0.41) 


0.65(0.33) 


2(2) 


-1.25(0.74) 


-5.41(0.61) 


-2.44(0.37) 


-7.78(0.39) 


O/I 




Y38 


22:19:02.26 


+56:06:36.5 


320 


3.7(1.7) 


0.48(0.01) 


2.28(0.10) 


6.25(0.26) 


1.02(0.34) 


3(2) 


-1.25(1.46) 


-5.58(0.97) 


-2.41(0.46) 


-8.06(0.43) 


O/I 




Y39 


22:19:02.62 


+56:06:17.3 


706 


4.2(1.7) 


0.47(0.01) 


2.26(0.11) 


6.13(0.33) 


0.44(0.25) 


2(5) 


-1.05(2.41) 


-5.23(2.22) 


-2.61(0.42) 


-7.79(1.71) 


O/I 




¥40 


22:19:02.77 


+56:05:40.6 


116 


14.9(1.2) 


0.48(0.01) 


5.03(0.53) 


6.27(0.12) 


1.85(0.27) 


20(10) 


-1.25(0.60) 


-5.66(0.48) 


-2.27(0.18) 


-8.16(0.31) 


O/I 




Y41 


22:19:03.80 


+56:05:10.3 


81 


12.1(1.2) 


0.47(0.01) 


2.27(0.06) 


5.85(0.36) 


0.68(0.22) 


4(2) 


-1.04(0.46) 


-5.26(0.28) 


-2.07(0.13) 


-7.21(0.24) 


O/I 




Y42 


22:19:04.36 


+56:05:50.8 


395 


7.1(1.2) 


0.48(0.01) 


2.29(0.08) 


5.94(0.41) 


0.74(0.36) 


11(12) 


0.00(0.86) 


-4.59(0.56) 


-2.04(0.48) 


-6.49(0.77) 


O/I 




¥43 


22:19:04.46 


+56:05:36.0 


1 


113.3 


0.48 


2.42 


6.09 


3.39 


36 


-2.98 


-8.68 


-2.20 


-7.59 


III 




Y44 


22:19:04.50 


+56:05:07.8 


215 


6.3(1.4) 


0.47(0.01) 


2.26(0.09) 


5.62(0.52) 


0.68(0.34) 


10(10) 


-0.32(0.70) 


-4.72(0.58) 


-2.19(0.43) 


-6.68(1.48) 


O/I 




Y45 


22:19:05.21 


+56:08:25.8 


10 


46.4(0.9) 


0.49 


2.31(0.03) 


5.42(0.04) 


0.85(0.08) 


7(1) 


-0.66(0.14) 


-4.92(0.14) 


-2.01(0.07) 


-7.22(0.08) 


O/I 




¥46 


22:19:05.72 


+56:04:52.7 


218 


6.7(1.3) 


0.48(0.01) 


2.27(0.09) 


5.30(0.59) 


2.07(0.69) 


62(37) 


0.18(0.60) 


-4.60(0.26) 


-1.42(0.37) 


-6.17(0.36) 


O/I 





Table 4 — Continued 



Name RA (J2000) DEC (J2000) Number 
(h m s) C ") 



log(d) 
log(kpc) 



A. 
(mag) 



log(Age) 
log(yr) 






L. 



log(M„,i.) log(M„v) log(Md,jt) log(Mrf,,,t) Stage 
log(Mo) log(Moyr-') log(Mo) log(Moyrl) 



Y47 22:19:05.79 +56:07:05.9 

Y48 22:19:05.82 +56:05:18.8 

Y49 22:19:05.94 +56:05:41.7 

¥50 22:19:06.02 +56:05:55.3 

Y51 22:19:07.65 +56:05:28.6 

Y52 22:19:08.03 +56:04:11.1 

¥53 22:19:08.30 +56:05:11.0 

Y54 22:19:08.45 +56:05:28.7 

Y55 22:19:09.38 +56:05:00.4 

¥56 22:19:09.39 +56:04:45.7 

Y57 22:19:09.70 +56:05:04.8 

Y58 22:19:09.84 +56:07:37.4 

¥59 22:19:09.94 +56:08:20.5 

Y60 22:19:10.94 +56:06:07.6 

Y61 22:19:12.07 +56:07:09.1 

¥62 22:19:13.80 +56:07:48.4 

Y63 22:19:15.06 +56:06:53.0 

Y64 22:19:15.18 +56:07:02.8 



138 


10.9(0.9) 


0.48(0.01) 


2.29(0.06) 


5.87(0.12) 


0.68(0.15) 


2 


-1.87(0.62) 


-6.15(0.31) 


-2.58(0.22) 


-8.21(0.31) 


III 


10 


108.1(0.6) 


0.47 


2.29(0.02) 


5.06(0.03) 


2.51(0.19) 


74(13) 


0.72(0.06) 


-4.08(0.04) 


-2.22(0.04) 


-7.37(0.08) 


O/I 


262 


4.9(1.6) 


0.48(0.01) 


2.27(0.09) 


5.89(0.31) 


0.71(0.31) 


5(8) 


-0.74(0.95) 


-4.83(1.21) 


-2.45(0.75) 


-7.32(2.15) 


O/I 


1 


41.0 


0.51 


2.00 


5.77 


3.58 


24 


-1.46 


-5.25 


-5.07 


-10.84 


O/I 


22 


14.5(1.6) 


0.48(0.01) 


2.23(0.06) 


5.68(0.30) 


3.09(0.40) 


67(25) 


-0.42(0.13) 


-4.60(0.15) 


-1.61(0.27) 


-6.47(0.42) 


O/I 


1 


130.6 


0.44 


2.00 


6.31 


4.80 


408 


-5.15 




-0.82 


-7.62 


III 


2 


387.4(0.6) 


0.45 


2.16(0.01) 


4.51(0.02) 


2.93(0.03) 


98(1) 


0.07(0.01) 


-4.41(0.01) 


-1.14(0.03) 


-5.57(0.03) 


O/I 


72 


9.4(1.8) 


0.48(0.01) 


2.27(0.07) 


5.51(0.35) 


1.28(0.43) 


18(10) 


0.27(0.49) 


-4.42(0.27) 


-1.61(0.18) 


-6.42(0.17) 


O/I 


1 


299.6 


0.46 


2.50 


4.12 


7.72 


607 


1.19 


-4.40 


-0.92 


-5.96 


O/I 


1 


193.3 


0.44 


2.00 


5.20 


7.24 


1003 


1.88 


-4.04 


-0.96 


-7.04 


O/I 


1 


30.7 


0.48 


2.00 


3.61 


5.91 


478 


1.22 


-2.92 


-2.32 


-6.35 


O/I 


1438 


4.1(1.5) 


0.48(0.01) 


2.27(0.11) 


5.72(0.54) 


0.36(0.24) 


2(3) 


-1.33(0.92) 


-5.42(0.53) 


-2.55(0.43) 


-7.76(0.88) 


O/I 


174 


8.4(1.1) 


0.48(0.01) 


2.27(0.07) 


5.45(0.40) 


0.32(0.14) 


2(2) 


-1.19(0.36) 


-5.28(0.25) 


-2.42(0.28) 


-7.33(0.52) 


O/I 


1134 


5.4(1.4) 


0.48(0.01) 


2.26(0.09) 


6.24(0.23) 


0.60(0.20) 


1 


-1.97(0.99) 


-6.00(0.68) 


-2.77(0.48) 


-8.48(0.47) 


III 


54 


9.8(1.2) 


0.47(0.01) 


2.20(0.07) 


5.76(0.18) 


0.66(0.14) 


3(2) 


-1.12(0.36) 


-5.31(0.25) 


-2.28(0.37) 


-7.15(0.45) 


O/I 


305 


15.9(1.2) 


0.49(0.01) 


2.82(0.42) 


6.75(0.04) 


2.18(0.12) 


26(7) 


-3.39(0.40) 


-8.72(2.28) 


-3.04(0.82) 


-10.24(0.57) 


III 


194 


4.0(1.7) 


0.48(0.01) 


2.24(0.10) 


5.92(0.47) 


0.49(0.24) 


1(1) 


-1.62(0.89) 


-5.73(0.75) 


-2.49(0.48) 


-8.20(0.44) 


O/I 


2239 


2.2(1.7) 


0.47(0.02) 


2.26(0.14) 


6.35(0.29) 


0.73(0.35) 


1(2) 


-1.85(13.16) 


-6.11(9.52) 


-2.85(0.85) 


-8.51(5.91) 


III 



I 



"Bad fit. 
''Evolved star? 



